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Abstract: Spontaneous parametric down conversion (SPDC) has been 
largely exploited as a tool for absolute calibration of photon-counting 
detectors, i.e detectors registering very small photon fluxes. In Q we 
derived a method for absolute calibration of analog detectors using SPDC 
emission at higher photon fluxes, where the beam is seen as a continuum 
by the detector Nevertheless intrinsic limitations appear when high-gain 
regime of SPDC is required to reach even larger photon fluxes. Here we 
show that stimulated parametric down conversion allow one to avoid this 
limitation, since stimulated photon fluxes are increased by the presence of 
the seed beam. 

© 2008 Optical Society of America 

OCIS codes: (270.4180) Multiphoton processes; (120.1880) Detection; (120.3940) Metrology. 



References and links 

1. H.A. Bachor, A guide to experimental quantum optics, (Wiley- VCH, New York, 1998). 

2. M. Genovese, "Research on hidden variable theories: A review on recent progresses," Phys. Rep. 413/6, 319-398 
(2005) and references therein. 

3. See, for a general review, N. Gisin, G. Ribordy, W. Tittel and H. Zbinden, "Quantum cryptography," Rev. Mod. 
Phys. 74, 145-195 (2000), and references therein. 

4. Giorgio Brida, Maria Chekhova, Marco Genovese, Alexander Penin, Ivano Ruo-Berchera, "The possibility of 
absolute calibration of analog detectors by using parametric down-conversion: a systematic study," J. Opt. Soc. 
Am. B 23, 2185-2193 (2006). 

5. B.Y. Zel'dovich and D.N. Klyshko, "Statistics of field in parametric luminescence," Sov. Phys. JETP Lett. 9 
40-44(1969). 

6. D.C. Bumham and D.L. Weinberg, "Observation of Simultaneity in Parametric Production of Optical Photon 
Pairs," Phys. Rev. Lett. 25, 84-87 (1970). 

7. D.N. Klyshko, "Use of two-photon light for absolute calibration of photoelectric detectors," Sov. J. Quant. Elect. 
10, 1112-1116 (1980). 

8. A. A. Malygin, A. N. Penin and A. V. Sergienko, "Absolute Calibration of the Sensitivity of Photodetectors 
Using a Two-Photon Field," Sov. Phys. JETP Lett. 33, 477-480 (1981). 

9. D.N. Klyshko and A.N. Penin, "The prospects of quantum photometry," Phys. Uspekhi, 30(8), 716-723 (1987). 

10. A. Migdall, "CoiTelated-photon metrology without absolute standards," Physics Today January, 41-46 (1999). 

11. G. Brida, M. Genovese and C. Novero, "An application of two photons entangled states to quantum metrology," 
J. Mod. Opt. 47 2099-2104 (2000). 

12. V. M. Ginzburg, N. G. KeratishviH, E. L. Korzhenevich, "Absolute meter of photodetector quantum efficiency 
based on the parametric down-conversion effect," G. V. Lunev, A. N. Penin, V. 1. Sapritsky, Opt. Eng. 32(11), 
2911-2916 (1993). 



13. G. Brida, M. V. Chekhova, M. Genovese. M. Gramegna, L. A. Krivitsky. and S. P. Kulik, "Conditioned Unitary 
Transformation on bipliotons," Pliys. Rev. A 70, 032332 (2004). 

14. G. Brida, M. Chekhova, M. Genovese, M. Gramegna, L. Krivitsky, M.L. Rastello, "Single-photon detectors cali- 
bration by means of conditional polarization rotation," J. Opt. Soc. Am. B 22, 488-492 (2005). 

15. G. Brida, M. Chekhova, M. Genovese. M. Gramegna, L. Krivitsky, M.L. Rastello, " Absolute quantum efficiency 
measurements by means of conditioned polarization rotation," IEEE Trans Instrum. Meas. 54 NO. 2, 898-900 
(2005). 

16. G. Brida et al., J. Mod. Opt. to appear 

17. K. Soda, I. Nishio, and A. Wada, "Analysis of noises in photocun'ent time-correlation spectroscopy of scattered 
light," J. Appl. Phys. 47, 729-735 (1976). 

18. R. G. Mclntyre, "Multiplication Noise in Uniform Avalanche Diodes", IEEE Trans. Electron. Devices, ED-13, 
164-168 (1966). 

19. The condition is fulfilled if the electronics is the same for the two detector and has a pass band much smaller than 
the band of the detectors. 

20. D.N. Klyshko, A. V. Masalov, "Photon noise: observation, squeezing, interpretation," Phys. Uspekhi, 38( 11), 
1203-1230 (1995). 



1. Introduction 

Prompted by the necessity of a precise calibration of photo-detectors both in the analog and 
photon-counting regime |[T]|2][3l, we presented recently |j4| a detailed theoretical analysis on 
the possibility to use the correlations of spontaneous parametric down conversion (SPDC) light 
ID |6l for calibrating analog detectors, as an extension of the technique developed in photon 
counting regime |I7][8l|9l[l0l[n][l2l[l3l[Ill[III- However, while the method is suitable for 
analog calibration at relatively low gain, which means a photon flux lower than lO"' photon/s, 
in higher-gain regime it is limited by the difficulty of collecting the same correlated modes in 
the two branches lH. 

In this paper we present a detailed theoretical analysis of a scheme based on the stimulated 
PDC (on the other hand a discussion of the uncertainty budget is left to a specifically addressed 
paper lfT6l ). We show that this scheme allow to overcome the problem mentioned above. Indeed, 
in this case the photon fluxes can be varied by varying the power of the coherent seed beam, 
without increasing the parametric gain G. On the other hand, the SPDC non-classical correlation 
at single-photon level which enables the absolute calibration of counting detectors survives, in 
some form, when a coherent seed is injected and the photon flux becomes macroscopic. 

2. Analog detection and quantum efficiency 

The schematic set-up for calibration of photodetectors by using PDC is shown in Fig. 1. We 
model the photodetection process in the analog regime as a random pulse train ifTTl 

n 

i.e. a very large number of discrete events at random times of occurrence f„. The pulse shape 
f{t) is determined by the transit time of charge carriers. We assume that /(f) is a fixed function 
with characteristic width Tp and a unit area. The pulse amplitude q„ is a random variable in 
order to account for a possible current gain by avalanche multiplication. The statistical nature 
of the multiplication process gives an additional contribution to the current fluctuations ifTSl . 
In an ideal instantaneous photocell response, without avalanche gain, all values q„ are equal to 
the charge e of a single electron and /(f) ^ 5(f). In the case of ideal quantum efficiency, since 
the probability density of observing a photon at time f at detector Dj (y = 1 , 2) is related to the 
quantum mean value {Fj{t)) of the photon flux, we calculate the average current output of Dj 
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Fig. 1. Scheme for absolute calibration of analog detectors by using stimulated PDC. The 
quantum efficiency of detector 2, collecting the seed beam and the stimulated emission, 
is estimated by the ratio between the cross-correlation function of the photocurrents and 
the auto-correlation of the photocurrent ii. 



as 

{'j)=LMit't„)} = f dt„{qj)fit~t„){Fj{t„)} 
n •' 

(1) 

where the factor {qj) is the average charge produced in a detection event. We have assumed the 
response function for the two detectors to be the same, f\ (f ) = fiit) = /(f)- 

At the same time, the auto-correlation and the cross-correlation functions for the currents can 
be expressed as 

{ij{t)ik{t + t)) ^ {qjnqkmf{t - tn)f{t - 1,„ + t)) 

= J J dtndt,n{qjqk)fit-tn)fit-tm + T){Fj{t„)Fk{tm)), 

(2) 

respectively for j = k where {Fj{tn)Fj(t,n)) is the auto-correlation function of the photon flux 
at detector j, and for j ^ k where {Fj{tn)Fk{tm)) is the cross-correlation between the fluxes 
incident on the two different detectors. By the substitution Fj = (Fj) + 5Fj, it is convenient to 
express them as 



{Fj{tn)mn)) = {Fj) {h) + {Fj)5[tn - t,n)5jk 

+ {■.dFj{t„)dFk{t,„):). (3) 

The second contribution, proportional to the photon flux when j = k, represents the the intrinsic 
and unavoidable component of the fluctuation that does not depend on the specific property of 
the field since it generates from the commutation relation of the quantum fields in the free 



space. The third one is the normally ordered correlation function of the fluctuation. Actually it 
is a function just of the difference t„ — t,,, and its typical variation scale provides the coherence 
time Tcoh of the PDC radiation. 

Now we introduce the quantum efficiency rjj of detector Dj, defined as the number of pulses 
generated per incident photon. In ID a real detector is modelled, as usual, with an ideal one 
(tj = 1) preceded by a beam splitter with the intensity transmission coefficient equal to the 
quantum efficiency of the real detector ll20l . Within this picture it is possible to take into account 
the quantum efficiency by the following substitutions: 

(Fjit)) ^ ri,{Fjit)) 
i-.Fjmit'):) iijrik{:Fjm{t'):). (4) 

Thus, (Fj{t)) being time independent, according to Eq. ([T]i we obtain: 

{ij) = rij{qj){Fj). (5) 



Eq. (|2|l becomes 



0(f)4(f + T)) = + r^j{qj)^ix){Fj)5jk 

fit-t„)fit-t„, + T)(: dFj{t„)dFtit„) :), 



(6) 



where we introduced the convolution of the response function of detectors ^(t) = 

Jdtfit)fit + T). 

3. Correlation functions of stimulated PDC 

Here we consider a type I PDC process. In the limit of monochromatic and plane-wave pump 
approximation, only pairs of modes with opposite transverse wave vectors, q and — q, and 
with conjugate frequencies, (Oi ~ COpump/'^- — and CO2 = cOpump/'^ + are coupled such that 
the energy and momentum conservation hold. The equations describing the down conversion 
process are the input-output transformations 



S?"'(q,n) = f/i(q,£2) S'l'iq,^!) +Vi{q,£l) ff,"'i-q~£L) 

-in 
'1 



fl|"'(q,n) = f/2(q,n) 4"(q,n) +y2(q,£2) 5;"^(-q, -i2) 



(7) 



linking the fields a^"' and d^"' after the non-linear interaction to the fields d'" and d'" before the 
interaction started. The coefficients f4(q,i2) and Vk{q,Q.) must satisfy the properties 

|f/,(q,n)|2-|y,(q,i2)|2 = l, (^=1,2) 
f/i(q,n)y2(-q,-i2) =C/2(-q,-n)yi(q,n), (8) 

which guarantee the conservation of the free-field commutation relations 

[5^-(q,n),5;«"t(q',n')] - d^.^d^.^^dnsi' 



[fl^"'(q,n),5^""(q',n')]-0. 

Here, Uk and V*. define the strength of the process and at the same time the range of trans- 
verse momentum and frequencies in which it takes place and thus are named gain functions. In 
fact, the finite length of the crystal introduces a partial relaxation of phase matching condition 
concerning the third component of the momenta. By selecting a certain frequency, the trans- 
verse momentum uncertainty, i.e. the angular dispersion of the emission direction, is not null 
{Aq ~ (/tantJ)^' for the non-collinear PDC, where I is the crystal length and ■& is the central 
angle of propagation with respect to the pump direction). On the contrary, by fixing the trans- 
verse momentum q, or equivalently a direction of propagation jJ, the spectral bandwidth AH is 
proportional to Z^'/^, and typically 1 /Ai2 ~ lO^^^s for type I. 

Hereinafter we will focus on the properties of the far field, observed at the focal plane of 
a thin lens of focal length /, placed at distance / from the crystal. The spatial distribution 
of the far field is, in this case, the Fourier transform of the field just at the output face of 
the crystal. This special imaging configuration is convenient to show the basic of calculation, 
nonetheless the validity of the results is more general. Thus, any transverse mode q is associated 
with a single point x in the detection (focal) plane according to the geometric transformation 
q — > 27rx/ [Xf). The far field operator in the space-temporal domain is therefore 

B,(x,f) SI- (^0x,i2) . (9) 

The mean value of the operator 4 (x,f) = Bj,(x,t)Bii{x,t) is 



on' 

("T'(0«.")«r(0x.n')). 

representing the intensity profile of the emission in the detection plane. This mean value has to 
be calculated over the following initial state, in which the field 1 is in the vacuum state whereas 
the field 2 is in a multi-mode coherent state: 

lr.«) = |O>i0|a(q,n))2, (11) 

where 

a(q,n) = a(q)5aoi2 

is the complex parameter associated to this state. The seed beam is therefore represented as a 
coherent state of field 2 with fixed frequency Q.() and a certain distribution of the transverse mo- 
mentum q. Here we assume that, for the modes that are stimulated by the seed, the spontaneous 
component of the emission is negligible with respect to the stimulated one. This corresponds to 
the assumption \ Vk{q,Q.o)\^ and, at the same time, |a(q)p ^ 1, that is a typical experimen- 
tal situation in which a few millimeters length crystal is pumped by a continuous pump and the 
seed has a power just around the microwatt or more. The intensities ( fTOb of the two stimulated 
beams after the crystal can be evaluated by using Eq.s (|7]), with the substitution q 2;rx/ {Xf). 
This leads to 



(/i(x,0)«|Vi(x,-no)|2|5(~x)|2 



(12) 



(/2(x,f)) ~ |C/2(x,no)|2|a(x)|2 

|V2(x,no)|2 + l)|a(x)|2 (13) 



where we defined 



Uk{x,a) = Uk[^x,Q. 

2n 
— s 
A/ 



yi(x,n) = yJ ^x,i2 



a(x) = a (^0x^ . (14) 

We used the first property of the gain functions in ([8j into the last line of Eq. ([T3j, in 
order to show explicitly the two contributions to the beam 2, one given by the original co- 
herent seed and the other from the stimulated emission, proportional to the parametric gain 
G = mflx|V]t(q,no)p. We notice that, aside from the spontaneous emission that is neglected, 
the generated beams conserve the original momentum distribution of the seed but weighted 
according to the gain function. 

The normal-ordered correlation function of the intensities fluctuation is defined as 



(:54(x,f)5/,(x',f'):) = (: 4(x,f)/,(x',0 :) - 

- (4(x,f))(/;-(x',f')) (15) 

The calculation leads to the following normal-ordered mean values: 

(:5/i(x,r)5/i(x',f'):)«0, (16) 

-i(Slo-Q.) (r-r')| 



(:5/2(x,0 5/2(x',f'):> ~ '|a(x 



n 



|y2(x,n)|2|[/2(x,no)|25. 



(17) 



Fluctuations of the stimulated field 1 are negligible, because they have the same order of mag- 
nitude of the spontaneous emission. Concerning the cross-correlation we have: 

(:5/i(x,f)5/2(x',f'):) « ^e-'("«+"H'-'')|5(_x)|2 

a 

f72(-x,no)C/|(-x,-£2) 

yi(x,-£2o)Vi*(x,n)5_,„/ 

+ C.C. (18) 

Let us consider two detectors, D\ and D2, to register the photons crossing the regions Ri and R2 
centered on the two correlated bright beams in the detection plane. Let the two areas be much 
larger than the dimension of the seed beam on the detection plane. This simply means to collect 
all the bright components of the emission, corresponding to the modes stimulated by the seed. 

We move now from the discrete modes representation of the electromagnetic field, in which 
it is artificially enclosed in a fictitious cube of side L, to a continuous set of modes for °o. 



Thus, in the following, the i2-sum and the Kronecker function 5^^/ in Eq.s ( fTTb and (fTSl ) are 
substituted respectively with the integral over D. and the Dirac function 5(x — x'). The fact that 
the second-order correlation function has an unphysical spatial behaviour, dominated by a delta 
function, comes from the starting assumption of a plane-wave pump, having infinite transverse 
dimension. In a realistic case the delta function is replaced by the Fourier transform of the pump 
transverse profile. 

The quantum mean values of the photon fluxes Fj = J]^Jj{x,t)dx reaching detectors 1 and 
2 can be obtained integrating equations ( fT2b and ( fTsT l respectively on Ri and R2- By applying 
relations ([8]) we have 

(Fi)^ f dx\Viix~ao)\^\ai~x)\\ (19) 

{F2)~J^ fl'x(|y2(x,no)|2 + l)|a(x)|2. (20) 

At the same time integrating Eq.s ( fTSI l and ( fTTI i over /?i x Ri and R2 x R2 respectively, we have 
the following auto-correlation functions for the photon fluxes 

(:5Fi(f)5Fi(f'):)«0 (21) 

{: 5F2it)5F2it') :) « J^^-mo-n) {,-,') 
a 

/ dx\V2{x,Q.)f\a{x)\^ + c.c., 

JR2 

(22) 

Here and in the following, we have assumed that in the limit of small G the gain functions 
\Uj{x,£l) \ ~ 1, according to the relations Then, integrating Eq. ( fTSI l over/?i x R2 gives the 
following cross-correlation: 

(: 8Fi{t)5F2{t') :) « / dx 

1 5(-x) (x, -£2o) yi*(x, + c.c. . 

(23) 

Although the seed beam is monochromatic with frequency £2o, the temporal width of both 
auto-correlation and cross-correlation in a fixed point x still depends on the spectral bandwidth 
AH of the gain functions. This explains why they originate from the correlation between the 
monochromatic seed and the spontaneous emission that is broadband. Therefore, the coher- 
ence time of the stimulated emission is still of the order of icoh = l/Ail, analogously to the 
spontaneous case. 

4. Quantum efficiency estimation 

The auto-correlation of the current fluctuations 5i\ (f ) = i\ (f ) — {i\ ) can be calculated by intro- 
ducing the result in Eq. ( |2T] ) in (|6]l, obtaining 



(5/i(05/i(f + T)) = 77i(^2^^(T)(Fi). 



(24) 



Under the condition of small gain (G <C 1) and large intensity of the seed (\a\^ ^ 1) the fluc- 
tuation of the current at detector Di is dominated by the shot noise component. 

By substituting Eq ( |23] | in under the same condition, for j = 1 and k = 2, one has 



(5/i(f)5/2(r + T))= 277i772 (?i)(^2)^(T)(Fi) (25) 

where we assumed Zp ^ Xcoh- We stress that it is the usual situation, the coherence time of 
SPDC being on the order of picoseconds or less and the typical resolving time of detectors on 
the order of nanosecond. In this case any fluctuations in the light power are averaged over Tp. 

The cross-correlation function of the fluctuations has the same form as the one obtained 
in the case of spontaneous down-conversion in our earlier work pl. This is the evidence that 
quantum correlations do not disappear when the emission is stimulated. In fact, the down- 
converted photons are still produced in pairs, although in the beam 2 the photons of the pairs 
are added to the bright original coherent beam propagating in the same direction. The factor 
2, appearing in Eq. ( |25] | can be interpreted as due to the fact that, for any down-converted 
photon of a pair propagating along direction 2, there is also the original photon of the seed 
that stimulated the generation of that pair. Formally it emerges from the time integration of two 
identical contribution in ( |23] |. the explicit one and its complex conjugated. 

Since we are interested in a relatively large power of incident light, we can at first consider 
detectors without internal gain, and assume that the charge produced in any detection event is 
equal to the single electron charge q, i.e. {q^) = q and {q^) — q^. Therefore, according to Eq.s 
(|24] | and (l25l l. the quantum efficiency can be evaluated as 

_l {dh{t)di2{t + T)) 

2{5k{t)5k{t + T)y ^^"^ 

For detectors in which electrons multiplication occurs the statistic fluctuations of this process do 
not allow to use Eq. ( l26l ) for absolute calibration. However, it can be performed by integrating 
Eq. i25[ over time T f?). It corresponds to evaluating the power spectrum of the fluctuations at 
frequencies around zero, namely much smaller than 1 /t,,. We would like to stress that in this 
case, the assumption fi (f ) — /2(f) — f{t) is not necessary. Since / dT,j^{T) = 1 we obtain 

I72W2 =X TTT , (27) 

2 ('i) 

In Eq. ( |27] | the statistic of the electron gain does not play any role. 



5. Conclusion 

With the purpose to eliminate some problems that appear when working with Spontaneous PDC 
we have extended our method to higher flux regimes, taking advantage of the stimulated 
PDC as a source of bright correlated beams. We show that this scheme allows one to overcome 
these former problem: an important result in view of metrology applications. Indeed this re- 
sult allows to link the photon counting regime with the photon rates of the traditional optical 
metrology supplying a single radiometric standard for all these regimes. Finally, we would like 
to stress once more that equations |26] and [27] between quantum efficiency and current fluctu- 
ations, even if identical to the ones obtained in our previous work Q, are here derived for a 
more general condition of stimulated emission. 



Acknowledgments 

The Turin group acknowledges the support of Regione Piemonte (ricerca scientifica applicata 
EI4) and San Paolo Foundation. The Russian group acknowledges the support of Russian Foun- 



dation for Basic Research (grant # 06-02-16393). Both acknowledge the joint project of Asso- 
ciazione Sviluppo del Piemonte by Grant RFBR-PIEDMONT # 07-02-91581-APS. 



